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We consider an anisotropic version of Baxter's model of 'sticky hard spheres', where a nonuniform 
adhesion is implemented by adding, to an isotropic surface attraction, an appropriate 'dipolar sticky' 
correction (positive or negative, depending on the mutual orientation of the molecules). The result- 
ing nonuniform adhesion varies continuously, in such a way that in each molecule one hemisphere is 
OO ■ 'stickier' than the other. 

>--^ ' We derive a complete analytic solution by extending a formalism [M.S. Wertheim, J. Chem. Phys. 

55, 4281 (1971) ] devised for dipolar hard spheres. Unlike Wertheim's solution which refers to the 
'mean spherical approximation', we employ a Percus-Yevick closure with orientational linearization, 
which is expected to be more reliable. 

We obtain analytic expressions for the orientation-dependent pair correlation function g{l,2). 
Only one equation for a parameter K has to be solved numerically. We also provide very accurate 
Cn ■ expressions which reproduce K as well as some parameters, Ai and A2, of the required Baxter factor 

correlation functions with a relative error smaller than 1%. We give a physical interpretation of the 

, , effects of the anisotropic adhesion on the p (1, 2). 

^1^ ■ The model could be useful for understanding structural ordering in complex fluids within a unified 

£t^ ' picture. 

'j^ ! PACS numbers: 61.20.Gy,61.20.Qg,61.25.Em 
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r^ , I. INTRODUCTION 

c. 

Anisotropy of molecular interactions plays an important role in many physical, chemical and biological processes. 
Attractive forces are responsible for the tendency toward particle association, while the directionality of the resulting 
bonds determines the geometry of the resulting clusters. Aggregation may thus lead to very different structures: in 
particular, chains, globular forms, and bi- or three-dimensional networks. Understanding the microscopic mechanisms 
underlying such phenomena is clearly very important both from a theoretical and a technological point of view. 
Polymerization of inorganic molecules, phase behaviour of non-spherical colloidal particles, building up of micelles, 

Zo , gelation, formation of a-helices from biomolecules, DNA-strands, and other ordered structures in living organisms, 

f^ ' protein folding and crystallization, self-assembly of nanoparticles into composite objects designed for new materials, 
are all subjects of considerable interest, belonging to the same class of systems with anisotropic interactions. 

f-^ • Modern studies on these complex systems strongly rely upon computer simulations, which have provided a number 

QQ ' of useful information about many properties of molecular fluids. 

f^ , Nevertheless, analytic models with explicit expressions for structural and thermodynamic properties still represent 

L| ' an irreplaceable tool, in view of their ability of capturing the essential features of the investigated physical systems. 

• ^H , At the lowest level in this hierarchy of minimal models on assembling particles lies the problem of the formation of 

^^ ■ linear aggregates, from dimers [l|,l4| up to polymer chains. This topic has been extensively investigated, through both 
J-j ] computer simulations and analytical methods. In the latter case a remarkable example is Wertheim's analytic solution 

- - I of the mean spherical approximation (MSA) integral equation for dipolar hard spheres (DHS), i.e. hard spheres (HS) 
with a point dipole at their centre [3| (hereafter referred to as I). For the DHS model, several studies predict chain 
formation, whereas little can be said about the existence of a fluid-fluid coexistence line, since computer simulations 
and mean field theories provide contradictory results 0, [a, [a, 0, I3- On the other hand, for mesoscopic fluids the 
importance of combining short-ranged anisotropic attractions and repulsions has been well established [9|, llO|, and 
hence the long-range of the dipolar interaction is less suited for the mesoscopic systems considered here, at variance 
with their atomistic counterpart. 

The aim of the present paper is to address both the above points, by studying a model with anisotropic surface 
adhesion that is amenable to an analytical solution, within an approximation which is expected to be valid at significant 
experimental regimes. 

In the isotropic case, the first model with 'surface adhesion' was introduced long time ago by Baxter [III, [131 ■ The 
interaction potential of these 'sticky hard spheres' (SHS) includes a HS repulsion plus a spherically symmetric attrac- 
tion, described by a square-well (SW) which becomes infinitely deep and narrow, according to a limiting procedure 
(Baxter's sticky limit) that keeps the second virial coefficient finite. 



14 , yj, [IJ, llll, M, lia, l20|, l2l| , 'sticky patches' [lO|, \M \M 
29j . The most common version of patchy sticky models 



Possible anisotropic variations include 'sticky points' [ic 
|24| . |25| . |26| . [23] and, more recently, 'Gaussian patches' [28 
refers to HS with one or more 'uniform circular patches', all of the same species. This kind of patch has a well-defined 
circular boundary on the particle surface, and is always attractive, with an 'uniform' strength of adhesion, which does 
not depend on the contact point within the patch |22| . 

In the present paper we consider a 'dipolar-like' SHS model, where the sum of a uniform surface adhesion (isotropic 
background) plus an appropriate dipolar sticky correction - which can be both positive or negative, depending on 
the orientations of the particles - yields a nonuniform adhesion. Although the adhesion varies continuously and no 
discontinuous boundary exists, the surface of each molecule may be regarded as formed by two hemispherical 'patches' 
(colored red and blue, respectively, in the online Figure 1). One of these hemispheres is 'stickier' than the other, and 
the entire molecular surface is adhesive, but its stickiness is nonuniform and varies in a dipolar fashion. By varying 
the dipolar contribution, the degree of anisotropy can be changed, in such a way that the total sticky potential can 
be continuously tuned from very strong attractive strength (twice the isotropic one) to vanishing adhesion (HS limit). 
The physical origin of this model may be manifold (non-uniform distribution of surface charges, or hydrophobic 
attraction, or other physical mechanisms), one simple realization being as due to an 'extremely screened' attraction. 
The presence of a solvent together with a dense ionic atmosphere could induce any electrostatic interaction to vanish 
close to the molecular surface, and - in the idealized sticky limit - to become truncated exactly at contact. 

For this model, we solve analytically the molecular Ornstein-Zernikc (OZ) integral equation, by using a truncated 
Percus-Yevick (FY) approximation, with orientational linearization (FY-OL), since it retains only the lowest order 
terms in the expansions of the correlation functions in angular basis functions. This already provides a clear indication 
of the effects of anisotropy on the adhesive adhesion. 

The idea of an anisotropic surface adhesion is not new. In a series of papers on hydrogen-bonded fluids such a 
water, Blum and co-workers |30l.l3ll.l32l| already studied models of spherical molecules with anisotropic pair potentials, 
including both electrostatic multipolar interactions and sticky adhesive terms of multipolar symmetry. Within appro- 
priate closures, these authors outlined the general features of the analytic solutions of the OZ equation by employing a 
very powerful formalism based upon expansions in rotational invariants. In particular, Blum, Cummings and Bratko 



[331 obtained an analytic solution within a mixed MSA/FY closure (extended to mixtures by Frotsykevich [3j|) for 
molecules which have surface adhesion of dipolar symmetry and at most dipole-dipole interactions. From the physical 
point of view, our model - with 'dipolar-like' adhesion resulting from the sum of an isotropic plus a dipolar term - 
is different and more specifically characterized with respect to the one of Rcf. [32], whose adhesion has a simpler, 
strictl y 'd ipolar', symmetry. From the mathematical point of view, however, the same formalism employed by Blum 
et al. [Sj] could also be applied to our model. Unfortunately, the solution given in Ref. [32] is not immediately usable 
for the actual computation of correlation functions, since the explicit determination of the parameters involved in 
their analytical expressions is lacking. 

In the present paper we adopt a simpler solution method, by extending the elegant approach devised by Wertheim 
for DHS within the MSA closure [3|, and, most importantly, we aim at providing a complete analytic solution - 
including the determination of all required parameters - within our FY-OL approximation. 

The paper is organized as follows. Section II defines the model. In Section III wc recall the molecular OZ integral 
equation and the basic formalism. In Section IV we present the analytic solution. Numerical exact results for some 
necessary parameters, as well as very accurate analytic approximations for them, will be shown in Section V. Some 
preliminary plots illustrating the effects of the anysotropic adhesion on the local structure are reported in Section VI. 
Fhase stability is breafly discussed in Section VII, while final remarks and conclusions are offered in Section VIII. 



II. HARD SPHERES WITH ADHESION OF DIPOLAR-LIKE SYMMETRY 

Let the symbol i = (r^, il^) (with « = 1, 2, 3, . . .) denote both the position r^ of the molecular centre and the orienta- 
tion fii of molecule i; for linear molecules, ^i = {9i, (pi) includes the usual polar and azimuthal angles. Translational 
invariance for uniform fluids allows to write the dependence of the pair correlation function g (1, 2) as 

(1,2) = (ri2,r!i,r!2) = (r,rJi,fi2,ri2) = (r, r!i, fia, ^i,), 

with ri2 = r2 — ri, r = |ri2|, and fir being the solid angle associated with ri2 = ri2/r. 

In the spirit of Baxter's isotropic counterpart [ll|, [S^, our model is defined by the Mayer function given by 

/SHS(i^ 2) = /HS(^) + t 6(1, 2) aS {r - a) , (2.1) 



where /^^(r) = (r — cr) — 1 is its HS counterpart, Q is the Heaviside step function {Q{x < 0) = 0, Q{x > 0) = 1) 
and 6 {r — a) the Dirac delta function, which ensures that the adhesive interaction occurs only at contact (cr being 
the hard sphere diameter). An appropriate limit of the following particular square well potential of width R — a 



$=^*(1,2) = 



can be shown to lead to Eq. (|2.ip . 

The angular dependence is buried in the angular factor 
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including the dipolar function 



e(l,2) = l+ai:'(l,2). 



£1(1, 2) = D{n-y,n2, nr) = 3(ui • f)(u2 • f) - m • U2 



(2.2) 



which stems from the dipole-dipole potential (/)'^'p~'^'p(1, 2) = —fj,^D{l, 2)/r^ {n is the magnitude of the dipole moment) 
and is multiplied by the tunable anisotropy parameter a. In the isotropic case, a = 0, one has e(l,2) = 1. Here and 
in the following, f coincides with f 12 = — i"2i , while u^ is the versor attached to molecule i (drawn as yellow arrow in 
Figure 1) which completely determines its orientation fli. Note the symmetry D{2, 1) = £'(1,2). 

The condition e(l,2) > must be enforced in order to preserve a correct definition of the sticky limit, ensuring 
that the total sticky interaction remains attractive for all orientations, and the range of variability —2 < D{1, 2) < 2 
yields the limitation < a < ^ on the anisotropy degree. The stickiness parameter t - equal to (12r)~ in Baxter's 
original notation [ll| - measures the strength of surface adhesion relatively to the thermal energy ksT (fc^ being the 
Boltzmann constant, T the absolute temperature) and increases with decreasing temperature. 

If we adopt an 'inter-molecular reference frame' (with both polar axis and cartesian z-axis taken along ri2), then 
the cartesian components of f and u^ are (0, 0, 1) and (sinfl^ coscp,;, sinOi sincp,;, cos^i), respectively, and thus 



D(l, 2) = 2 cos 01 cos 6*2 — sin^i sin02 cos {ipi — (^2) 



(2.3) 



The strength of adhesion between two particles 1 and 2 at contact depends - in a continuous way - on the relative 
orientation of Ui and U2 as well as on the versor ri2 of the intermolecular distance. We shall call parallel any 
configuration with Ui • U2 = 1, while antiparallel configurations are those with Ui • U2 = — 1 (see Figure 1). For all 
configurations with D{1, 2) > 0, the anisotropic part of adhesion is attractive and adds to the isotropic one. Thus, the 
surface adhesion is maximum, and larger than in the isotropic case, when Ui = U2 = ri2 and thus e(l, 2) = l+2a (head- 
to-tail parallel configuration, shown in Figure lb). On the contrary, when Z?(l,2) < the anisotropic contribution 
is repulsive and subtracts from the isotropic one, so that the total sticky interaction still remains attractive. Then, 
the stickiness is minimum, and may even vanish for a = 1/2, when Ui = — U2 = ri2 and thus e(l,2) = 1 — 2a 
(head-to- head or tail-to-tail antiparallel configurations, reported in Figure Ic). The intermediate case of orthogonal 
configuration (u2 perpendicolar to Ui) corresponds to D(l, 2) = 0, which is equivalent to the isotropic SHS interaction. 

It proves convenient to 'split' /^^^(1,2) as 



/^"^(l,2) = /o(r)+/e.(l,2). 



(2.4) 



/(,(r) = /HS(r) +ta5{r-a) = /-°shS(^) 
/ex(l,2) = (ai) a5{r~a) D{1,2), 



(2.5) 



where the spherically symmetric fo{r) corresponds to the 'reference' system with isotropic background adhesion, while 
/ox(l,2) is the orientation-dependent 'excess' term. 

We remark that, as shown in Ref. I (see also Table I in Appendix A of the present paper), convolutions of Z^^^- 
functions generate correlation functions with a more complex angular dependence. Therefore, in addition to D{1, 2), 
it is necessary to consider also 



A(l, 2) = Ui • U2 = cos 01 cos 02 + sin 01 sin 02 cos (ipi — if2) , 



(2.6) 



where the last equality holds true in the inter-molecular frame. The limits of variation for A(l, 2) are clearly —1 < 
A(1,2)<1. 



III. BASIC FORMALISM 

This section, complemented by Appendix A, presents the main steps of Wertheim's formalism, as well as its 
extension to our model. 

A. Molecular Ornstein-Zernike equation 

The molecular OZ integral equation for a pure and homogeneous fluid of molecules interacting via non-spherical 
pair potentials is 

h{l, 2) = c(l, 2)+p I dr, { c(l, 3) /i(3, 2) )^,^ , (3.1) 



where h{l,2) and c(l,2) are the total and direct correlation functions, respectively, p is the number density, and 
f/(l,2) = 1 + h{l,2) is the pair distribution function [3y, [33, [3^. Moreover, the angular brackets with subscript fi 
denote an average over the orientations, i.e. {■ ■ ■)fi = (47r) J dfl ■ ■ ■ . 

The presence of convolution makes convenient to Fourier transform (FT) this equation, by integrating with respect 
to the space variable r alone according to 

F (k, fli, fla) ^ [ dr F{r, ^i, f^a) exp(ik • r). (3.2) 

The r-space convolution becomes a product in k-space, thus leading to 

h{k,ni,n2)=dik,ni,n2) + p (dik,ni,n3)h{k,n3,n2)) . (3.3) 



As usual the OZ equation involves two unknown functions, h and c, and can be solved only after adding a closure, 
that is a second (approximate) relationship among c, h and the potential. 

B. Splitting of the OZ equation: reference and excess part 



The particular form of our potential, as defined by the Mayer function of Eq. (|2.ip . gives rise to a remarkable exact 
splitting of the original OZ equation. Using diagrammatic methods [36l. 1371. l38l| it is easy to see that both c and h can 
be expressed as graphical scries containing the Mayer function / as bond function. If /^^^ = /o + /ex is substituted 
into all graphs of the above series, each diagram with n /-bonds will generate 2" new graphs. In the cluster expansion 
of c, the sum of all graphs having only /o-bonds will yield co{r) = c"'°^^^(r), i.e. the the direct correlation function 
(DCF) of the reference fluid with isotropic adhesion. On the other hand, all remaining diagrams have at least one 
/ex-bond, whose expression is given by Eq. (|2.5p . Thus, in the sum of this second subset of graphs it is possible to 
factorize at, and we can write 

cSHS(i,2) = co(r)+Cex(l,2), (3.4) 

Cex(l,2) = (at) ct(l,2). ^^-^^ 



Similarly, for h we get 



/iSHS(i^2) = /io(r)+/iex(l,2), (3.6) 

/iex(l,2) = (at) /^t(i,2). ^■'■') 



Note that this useful separation into reference and excess part may also be extended to other correlation functions, 
such as 7(1, 2) = h{l,2)-c{l,2),g{l,2) = l + h{l,2), and the 'cavity' function y(l, 2) = g(l, 2)/e(l, 2). The function 
7 coincides with the OZ convolution integral, without singular (5-terins. Similarly y is also 'regular', and its exact 
expression readsy (L 2) = exp [ 7 (1, 2) + B{1, 2) ], where the 'bridge' function B is defined by a complicated cluster 
expansion [36|, l37|, |38[ . 

From Eqs. p.4p - p.7p . which are merely a consequence of the particular form of /ex in the splitting of /®^^, one 
immediately sees that, if the anisotropy degree a tends to zero, then 

lim ccx(l, 2) = lim /iex(l, 2) = lim yex(l, 2) = 0. (3.8) 

a — >0 a — >0 a — >0 

Note that the spherically symmetric parts cq and ho must be related through the OZ equation for the reference 
fluid with isotropic adhesion {reference OZ equation) 

ho{r) ^ co(r) + p I dr^ coir^) /io(?'32)- (3.9) 



Thus, substituting c and h of Eq. (j3.ip with cq + Cox and /iq + /lox, respectively, and subtracting Eq. (|3.9|) . wc find 
that Ccx and ft-cx must obey the following relation 

/iex(l,2) = Ccx(l,2) + p /dra [ coim) {h^43,2))^^ 

+ { Cex(l, 3) ),j3 ho{r32) + { Ccx(l, 3) /Jex(3, 2) )^^ ] . 

and when 

(Cex(l,3)),,, =(/iex(3,2)),,^=0 (3.10) 

the orientation-dependent excess parts Cex and ft,ox satisfy the equality 

/lex(l,2)=Cex(l,2)+p /dra (Cex(l,3)/lex(3,2))^ , (3.11) 



which is decoupled from that of the reference fluid and may be regarded as an OZ equation for the excess part {excess 
OZ equation). As we shall see, condition (j3.10p is satisfied in our scheme. 

We stress that, in principle, the closures for Eq. (|3.9[) and Eq. (|3.1ip . respectively, might be different. In addition, 
although the two OZ equations are decoupled, a suitably selected closure might establish a relationship between Fq 
andF {F = c,h). 



Percus-Yevick closure vi^ith orientational linearization 



For hard-core fluids, h and c inside the core are given by 

h{l,2) = -1 for 0< r < cr. 



c(l,2) = -[l + 7(l,2)] for 0<r<cr. ^^-^^^ 

At the same time, we have the following exact relations 

h{l,2) = ,9(1,2)- 1 = 6(1,2)2/(1,2) -1, 

c(l,2) = /(l,2)[l+7(l,2)]+e(l,2)[y(l,2)-l-7(l,2)]. 

Since c, h and g are discontinuous for hard-core fluids and involve (5-terms for sticky particles, it is more convenient 
to define closures in terms of y and 7, which are still continuous and without (5-singularities. The Percus-Yevick 
approximation for molecular fluids with orientation-dependent interactions corresponds to assuming 

/^ (1, 2) = 1 + 7(1, 2) everywhere, (3.13) 

and thus, for the DCF, 

cPY(l,2) = /(l,2) [1 + 7(1,2)], (3.14) 



which impUes that c vanishes beyond the range of the potential. 

However, the dependence of 7(1,2) on angles may still be very complex. A possible procedure is to perform a 
series expansion of all correlation functions in terms of an infinite set of rotational invariants, which are angular basis 
functions - related to the spherical harmonics ~ having the property of rotational invariance valid for homogeneous 
fluids [33| . Unfortunately, the full PY approximation requires an infinite number of expansion coefficients for both 
c(l, 2) and /i(l, 2). This approach is usually impracticable, but sometimes even unnecessary, as it is possible that the 
most significant angular basis functions arc included in a small closed subset of that infinite set. Indeed this happens, 
for instance, in the DHS model within the MSA [3|, where the set {1,A,D} is the required subset. Although this 
does not happen in our model, we shall argue that the same truncation is sufficient due to the dipolar symmetry of 
the anisotropic adhesion. 

Indeed, a natural assumption is that the only nonzero harmonics in c(l, 2) and /i(l, 2) are those contained in /(I, 2) 
and those which can be obtained from that set by convolution [30|. Now, the angular basis functions included in our 
/-bond are only 1 and £), but the convolution of two /-bonds involves the angular average of two D's, which yields 

1 



-[D{k,ni,n2) + 2A(k,ni,n2) 



in k-space, and thus generates also A. Consequently, we will expand any angle-dependent correlation function F as 

F{1,2) = Foir) + FA{r)A{l,2) + FD{r)D{l,2) + ■ ■ ■ , (3.15) 

neglecting all higher order terms. In other words, we assume that all angular series expansions can be truncated after 
these first three terms, linear with respect to the angular basis functions. Using this spirit in the PY approximation, 
given by Eq. p.l4p . we obtain the following PY correlation functions with orientational linearization (OL): 



PY-OL 



(1,2) ^ co(r)+CA(r)A(l,2)+ci3(r)i?(l,2) 

= coir) + {at) [4(r)A(l,2)+4(r)i?(l,2) 



(3.16) 



and 



PY-OL 



(1,2) = ho{r) + h^{r)A{l,2) + hDir)D{\,2) 

= ho{r) + iat) \hi{r)A{l,2)+h^^{r)D{l,2) 



where 



co(r) = Ao ad {r - a) 
ca(»') = Aa cr6{r- a) 
CD{r) = Ad cr5{r - a) 



for r > a, 



(3.17) 



(3.18) 



with 



Ao = yn^r) t 
Aa ^yr{<j)t 
Ac = [ y£^(CT) + a y^^{<j)] t = yl^{a) t + aAo, 



(3.19) 



j/r(r) = l + 7o(r) 

yA^ir) = 7A(r) == (at) yl^ir) 

yo^ir) =lD{r) = (ai) yoi^)- 

Clearly for /(I, 2) no truncation is required, as the expansion 

/o(r) = /'-SHS(^) ^ ^HS(^) +ta5{r~a) 

/A(r) = 

foir) = [at) a5{r - a) 



(3.20) 



(3.21) 



is exact. 



It can be shown that c(l,2) and /i(l,2) must have the same approximate form in view of the OZ equation. Eq. 

(EH). 

The solution of the original OZ equation (|3.ip is then equivalent to the calculation of the radial coefficients 
co('')i ca(?'), C£)(r) and ho{r),h/^{r),h]:)(r), which are the projections of c(l,2) and /i(l,2) onto the angular basis 
{1, A, £)}. The core condition on h, Eq. (|3.12p . becomes 

hoir) = -I 

h\{r) = ). for < r < cr. (3.22) 

hD{r) = 

Note that in the zero density limit 7(1, 2) = p J dr^ ( c(l, 3) h{3, 2) )^ must vanish, and thus y^^{l, 2) -^ 1, i.e. 

\hnyr{r) = 1, Hm J/^M = lim y^M = 0, 

p— >0 p— >0 p^O 

while both c (1, 2) and /i (1, 2) must reduce to /(I, 2): 

limp_o-Fo(r-) = /oC?-) 
limp^o^AW=0 (i^ = c,/i), 

hmp^nFoir) = /^(r) 



and 



lim Ao = t, lim Aa = 0, lim Ad = at. (3.23) 

p^O p^O p^O 

Moreover as a ^ all A- and D-coefficients of c, h and y vanish, so that the isotropic adhesion case is recovered. 
Finally, it is also worth stressing that the same J-tcrm arises in c, h and g, that is 

F(l,2)-F,eg(l,2) + f,ing(l,2) iF = c,h,g), 

where Frcg is the 'regular' part (i.e., the part with no ^-singularity, and - at most - some step discontinuities), while 
Fsi„g(l,2) = aS {r — cr) A(l,2) is the singular term representing the anisotropic surface adhesion ( with A(l,2) = 
Ao + AaA(1,2)+Az3Z?(1,2) ). 

D. Integral equations for the projections of c and h 

In the following, we extend Wertheim theory [3| to our model, in order to obtain the radial projections of c and h. 
The PY-OL approximation to the excess anisotropic part of the correlation functions is 



cP/-OL(i, 2) = CA(r) A(l, 2) + cz,(r)i?(l, 2) 
CY-OL(i, 2) = hA{r)A{l, 2) + hn{r)D{l, 2), 



PY-OLn r,\ _ u . f^\/\f^ o\ 1 h^f^\Tifi o^ (3-24) 



thus verifying the required property ( Cox(l,3) )q = ( /iox(3, 2) )q = described in Section III, and allowing the 
splitting of the molecular OZ equation into a reference and an excess part. 

The fir st p art is the reference PY equation, and coincides with that solved by Baxter for the fluid with isotropic 
adhesion [l3,[i3|: 

/io(r) =Co(r) + p (ho-kco) 

hair) = -1 < r < cr (3.25) 

Co(r) = Ao cr (5 (r — ct) r > a, 

where the symbol • denotes spatial convolution, i.e. {A-k B){ri2) ~ J A{ri3)B{r32) dr^. 

The second part is the excess PY-OL equation, given by Eq. (|3.1ip coupled with the PY-OL closure. Following an 
extension of Wertheim's approach, as described in detail in Appendix A, Eq. (|3.1ip can be splitted into the following 
system for the A— and D— projections of c and h : 

h^ir) = CAir) + ^p ( CA*/iA + 2 c2)*/ic ) .„ ^g) 



hoir) = c%{r) + Ip { ca* h% + c%* Ha + c\ 



D 



where c^jj{r) and h^{r) are defined by the relationship 

F%{T)^FD{r)~i r^^^dx {F = c,h), (3.27) 

Jr ■^ 

whose inverse is Q 

Foir) - F%{r) -\ f F^ (x) x^ dx. (3.28) 

*" Jo 



The eore conditions become 



titr) = -3K } f°^- 0<'^<'^' 



(3.29) 



with 



oo 



if = / '^^^^ dx = ifrcg + Ai5 , (3.30) 



K,,^ = / ^^^^£iW ^^_ (3 3^) 



Note that the presence of the (5-singularity in ho{x) requires the specification of a as lower integration limit, unlike 
the case of Ref. I where only the regular part Krcg is present. Moreover, since /i £>(?') = at h]j{r), from Eq. ()3.30p 
one could also write 

K = atlC , (3.32) 

which shows that K is related to the anisotropy degree, and vanishes both in the symmetric adhesion case {a = 0) 
and in the HS limit (t = 0). Since ho{r) — > foif) = (ctt) ad{r — a) in the zero density limit, one then finds that 

lim K = at. (3.33) 

Finally, the PY-OL closure for the new DCFs reads 

c^ir)=A^aS{r-a)\ 

d[j(r) = Ad (JO (r - a) j "■ ' 

(for simplicity, here and in the following we omit the superscript PY-OL). 

E. Decoupling of the integral equations 

It is possible to decouple the two equations for A- and ZJ-coefficients, by introducing two new unknown functions, 
which are linear combinations of the previous ones. As shown in Appendix A, if we define Fi{r) and F2{r) [F = c, h) 
through the relations 

Fi(r) = (3A)-^[FA(r)-F0(r)] .p ^ . ^ 

F,{r)^i3C,)-'[FMr)+2FUr)] ^ "' ' ' 

then we get the OZ equations 

hi{r) = ci(r) + pi (hi^ci) 
h2{r) = C2(r) + p2 (/i2 * C2), 

with the following densities and core conditions 



Pi 

P2 



Ci p 
C-2 P, 



h2{r) = -2K/C2 



for < r < a. 



The decoupling of the three different projections of c and h is remarkable: the molecular anisotropic OZ equation 
reduces to a set of three radial integral relations, which may be regarded as OZ equations for three 'hypothetical' 
fluids (labelled as 0, 1, 2) with spherically symmetric interactions. We stress that there is not a unique solution to the 
decoupling problem, since - in principle - there exist infinite possible choices for (£i,£2)- The final results are clearly 
independent of the values of (£i,£2)- 

In the present paper, we adopt Wertheim's choice, i.e. £1 = —K and £2 = 2/\ , which leads to 



F2{r) = S[FUr) + ^F^{r) 



[F^c.h), 



(3.35) 



Pi = -Kp 
P2 = 2Kp, 



hi{r) = -l 
h2{r) = -1 



for < 7' < (7 



(3.36) 



(in Ref. I, Fi and F2 were denoted as F- and F+, respectively). 

Note that the auxiliary fluids have densities different from that of the reference fluid (the negative sign of pi poses 
no special difficulty). 

We can also write 



with 



Fm (r) = Fm,^cg{r) + A™ aS (r - ct) , 



(3.37) 



^l,reg(r)=3i^[F^,,,g(r)-FA,.eg(r)] 
i^2,reg (r) = 3^ [F^,„g(0 + I^A.reg(r)] 



and 



Ai = 3!^ (Ad - Aa 



= 3K [^-D,rcg(o-+) - /lA,rog(a-'*')] t + a g^^ Aq 

A2 = ^ {^D + ^Aa) 

= JK [^D,rcg{<y^) + 5/lA,rcg(CT+)] t + a ^ Kq. 

(since ^...{<y) = /i...,rcg(CT^) — c...^i.cg(CT^), and c...,rcg(o'+) = within the PY-OL closure). 

Knowing the correlation functions Fi{r) and F2(r) (with F = c,h), one can derive Fa (r) , F'q (r), i.e. 

Fa (r) = 2K [F2 (r) - Fi (r)] 
F" (r) = 2K [F2 (r) + iFi (r)] , 



and 



Aa = 2A'(A2-Ai) 
Ad = K{2A2 + ^i) ■ 

Finally, from Fa (r) , F^ (r) one has to evaluate Fa (r) , Fq (r), by employing Eq. p.28p . We note the 
points: 

i) Insertion of h^{r) — /i^ i-cA^) + Ad cr6{r — a) into Eq. p.28p yields hoir) = /iD,reg(?') + A^ (Td{r — a] 



(3.38) 



(3.39) 



At r = 2ct /iDaog and /i^j r^g have the same discontinuity. We also get 



r > (7. 



(3.40) 



(3.41) 

following 
, with 

(3.42) 
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/lD,reg(a+) = /l2,,,.eg(T+) + 3K,eg. (3.43) 

Clearly, these results must agree with those obtained from Eq. p.27p . i.e. 

ipir) = J^ hoix) X ^ dx ^ Ad {(7 - r) + f^ hD.icgix) x ^ dx. 

In order to recover Eq. (|3.43p along this second route, note that '0(r) is not continuous at r = cr. In fact, from Eqs. 
([3:30)) and ((33T|1 follows ^(cr") = K whereas 'il^{a+) ^ A'„g. 

ii) Similarly, for cjjlr) wc obtain c^lr) ~ CD,-cg{r) + A^ ad{r — cr), with 



CDscg{r) =c23 (r) -3r' 



CD.rciri^)^^ dx + ^D(7^ 0{r - a) 



(3.44) 



since /„ ^{x — a)x'^dx — a^9{r — a). On the other hand, from Eq. p.27p one easily finds that 

CD{r) =c%{r) for r > a. (3.45) 

iii) By applying the relationship p.28p to cd (r) , using Eq. (|3.45p and noticing that cq {x) = for r > ct within the 
PY-OL approximation, leads to a sum rule: 

/ c%{x)x'^dx= / c%^^g{x)x'^ dx + Aoa^ =0, (3.46) 

Jo Jo 

that we will exploit later. 

IV. ANALYTIC SOLUTION 

We have seen that the molecular PY-OL integral equation (IE) for our anisotropicSHS model splits into three lE's 

h^{r) = -l 0<r<(T (m = 0,l,2), (4.1) 

Crn{r) = A„i aS{r - a) r > a 

where 

Pi = -Kp (4.2) 

, P2 = 2Kp, 

and the 'amplitudes' of the adhesive (5-tcrms are 



Ao = [1 + /io,rcg(CT+)] t = yl^{a)t 

Arn - /im.rcg(fT+) t +V= V^ {a)t + V (m = 1, 2) , 



(4.3) 



with 



1 at vl^ja) Xrcg ^ ^ 1 « Aq i^ ^ .^ ^^^ 

3 A' if 3 if A' ■ ^ ' ' 

Here, the new expressions of Ai and A2 have been obtained from Eqs. (|3.39p with the help of Eqs. (|3.43p and (|3.38p . 
The essential difference with respect to Ref. I lies in the closure, which is - of course - related to the model 
potential. While Wertheim's paper on DHS [31 employed the MSA closure, which performs properly for long-ranged 
electrostatic potentials at low strength of interaction, our PY-OL closure is more appropriate for the short-ranged 
potential of the present model. 
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The first integral equation lEO is fully independent, whereas lEl and IE2 depend on the solution of lEO (unlike 
the case of Ref. I) , because of the presence of Ag inside Ai and A2 . While lEO is exactly the PY equation for the 
reference SHS with isotropic adhesion solved by Baxter [111, [l3l, lEl and IE2 are different from both Wertheim's 
MSA solution for DHS and Baxter's PY solution for SHS. We remark that the closures for lEl and IE2 are not PY 
as Ai and A2 - given by Eq. (|4.3[) - differ, by the term V, from those appropriate for the PY choice, corresponding 

toAr=yrwi- 

Consequently, lEl and IE2 can be rekoned as belonging to a class of generalized PY (GPY) approximations, 
introduced in Ref. [39], which admit an analytic solution. Thus, the PY-OL closure for c(l, 2) leads to a PY integral 
equation for co{r), coupled to a two GPY integral equations for ci(r) and €2(7") (which are linear combinations of 
CA(r) and c%{r)). 

On comparing the three lE's and their closures given by Eq. (|4.ip . it is apparent that they have exactly the same 
form, but differ by the density pm and the expression for A^. The first integral equation lEO corresponds to an 
isotropic SHS fluid with density p. On the other hand, lEl and IE2 refer to 'auxiliary' isotropic SHS fluids with 
densities pi and p2, and adhesion parameters Ai and A2, respectively. Note that, according to Eqs. (|3.39p , Am is not 
evaluated at the actual density pm of the auxiliary fluid, but at the real density p . These remarks strongly suggests 
that the solutions of lEO, lEl and IE2 can be expressed in terms of an unique solution - the PY one for isotropic SHS 
- by changing only pm and Am- This can be achieved by the formal mapping 

Fiir) = F'^°^'^^{r;r]i,Ai) {F = q,c,h) (4.5) 

F2(r)=F'-SHS(^.^2^A2), 

where 770 = 77 is the real volume fraction, while 771 and 772 are 'modified volume fractions' of the 'auxiliary' fluids 1 
and 2, i.e., 

770 = 7/ = (tt/G) pa^ 

771 = -Kt] (4.6) 

772 = 2K7] . 

In Eqs. (|4.5p q{r) denotes the Baxter factor correlation function, introduced in the next Subsection. 

It is worth noting that this result for SHS mirrors the analog of the MSA solution for DHS [3| where all the three 
harmonic coefficients of can be expressed similarly, in terms of a single PY solution for the reference HS fluid. 



A. Baxter factorization 

We shall now solve Eqs. (|4.ip by using the Wiener-Hopf factorization due to Baxter [i2|. Let us recall its basic 
steps. After Fourier transforming the OZ equation for a one-component fluid with spherically symmetric interactions, 
one assumes the following factorization: 

l~pc{k)^Q{k)Q{~k), 
Q{k) ^ 1 - 2npJ^ q{r) e'^'-dr. ^ '> 

Then it can be shown that the introduction of the 'factor correlation function' q{r) allows the OZ equation to be cast 
into the form [l2| 

re (r) = -q'(r) + 27rp j^jiu q (u) q' {r + u) , ,^ ^. 

rh{r) ~ —q'{r) + 2iTp L du q{u) {r ^ u) h {\r — u\) , 

where the prime denotes differentiation with respect to r. Solving these Baxter equations is tantamount to determining 
- within a chosen closure - the function q{r), from which c{r) and h {r) can be easily calculated. It is also necessary 
to remember that, for all closures leading to c (r) = for r > a, one finds g (r) = for r > a [39|. 
On applying Baxter's factorization to Eqs. (|4.ip . we get 

rhm {r) = -q'mif) + 27rp„ / du q,n (u) (r - u) h^ {\r - u\) . (4.9) 



with m = 0, 1, 2. Now the closure Cmir) = A„i aS (r — a) for r > a implies that the same (5-tcrm must appear in 
hm (r)- Thus, for < r < a, using hm{r) = — 1 + Am crS (r — cr), we find 
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with 



9m (^) = «'"?' + b„icr - A,„ (7^(5 {r- a) ^ 



am= 1 - 27rp,„ /g'^ du q^^ (u) 
brna = 27rp,„ /g'^ dw (?,„ (w) u . 



(4.10) 



The (5-terni of g^„(r) means that Qml'') has a discontmuity qmio'^) — Qm{o'~) = — AmCr^, with qrn{o'~^) ~ 0. Integrating 
q'^(r), substituting this result into Eqs. (j4.10p . and solving the corresponding algebraic system, we find the following 
solution 



qm{r) 



\a^(r - uY + (am + &m) ^(r - ct) + A,; 



< r < cr, 
otherwise, 



(4.11) 



HS/ \ 1277m A^ 



1 -?7r) 
1 - ?7m 



6m = 6"^(?7m) + 

rim = (7r/6) pmO-^ 
a"^(x) 



(1 — x) 2(1 — x) 



(4.12) 

(4.13) 
(4.14) 

(4.15) 



From the first of Eqs. 
for < r < cr 



we get the DCFs Cra{r) = c,„_iog('') + A™ a5{r — a), where Cm.-[cg{T) = for r > a, and 



Cm,rcg{r) = --^Vm O^ f " j + 6?7m [ (ftm + ^m)^ - 2am Am ] (-] 



-flm - 12?7m A,2 



(4.16) 



The second of Eqs. (|4.8p yields the total correlation functions h„i{r) ~ /im,rcg(^) + A™ (T(5(7' — a). For ?- > cr, Eqs. 
(|4.9p becomes 



/J' " du qm (u) i?m,rcg {r - u) 

iJ„ ,„ (r) == 1277„ cr-^ i + ^r-a ^" 9m (u) {u - r) + Amcr^ ^m (r - (t) 

J^ du q„i (u) H,n^rcg {r-u) r > 2cr, 



a <r <2a, 



(4.17) 



where Hm (r) = rhm (r). Due to the last term of Eq. (j4.17p and the discontinuity of qm (r) at r — a, /im,rog(^) has a 
jmirp of at r = 2cr [40,|4l|: /im,rcg(2cr+) - /im,,og(2cr") = -6?7m A^. 



B. An important relationship 

In Appendix B it is shown that a remarkable consequence of the sum rule (|3.46p is the condition 

a2 == ai , (4.18) 

that will play a significant role in the determination of the unknown parameters Ai, A2 and K (see Appendix B). 
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C. Reference fluid coeflicients 

The 171 = case corresponds to Baxter's PY results for the reference fluid of isotropic SHS particles [lil, [i4|. We 
have: qo{r) = g'''°^^^(r;77, Aq), and 



Co(r) == c''°^^^r;T],Ao) = 4lf^^r;7^,AQ)+Ao aS{r - a) 
ho{r) = /i-°SHS(^.^^^^) ^ hf,f''^{r;i^,Ao) + Ao aS{r - a) 

(for simplicity, we omit ~ here and in the following - the superscript P Y) . 



D. A— and D— coeflicients 

We write Qmir) = q^^°^^^ {r; ijm, Am) with m = 1,2. Then, 

i) For the A-coejficients, after recalling Eq. (|3.37p and exploiting Eqs. (|3.4ip . wc end up with: 

CA(r) = 2A'[coa-cg(?-;2i4'77,A2)-co,rcg(r;-ii:?7,Ai)] + Aa (T6{r - a) 
hAir) = 2K [/io,rcg(?-; 2Kt^, As) - /loa-cglr; ~Kt^, Ai)] + Aa cr(5(r - cr). 

ii) For the D- coefficients ^ we get 



(4.19) 



(4.20) 



r c2,(r) = 2is:[co,reg(r;2i^?7,A2) + ico,rog(r;-ii:?7,Ai)] + Ad cr5{r - a) , . 

\ hlir) = 2K [/io,reg(r; 2K7j, A^) + Iho^rcgir; -Krj, Ai)] + A^ aS{r -a). ^ ' ' 

Finally, from c'jj{r) and h'jj{r) we can calculate coir) and hnir), as described by Eqs. (j3.44p and p.42p . respectively. 
In short, a) our PY-OL solution - {cq,ca,cd} and {ho,hA,h]:)} - satisfies both the PY closures and the core 
conditions; b) all coefficients contain a surface adhesive 5— term; c) {/iqi^Ai^d} a-U exhibit a step discontinuity at 
r = 2a. 



V. EVALUATION OF THE PARAMETERS K, Ai AND Aa 

The calculation of the Baxter functions g,„s (m = 0, 1, 2) requires the evaluation of A', Ai, and A2, for a given set 
of a, ?7 and t values, a task that we address next. 



A. Exact expressions 

Four equations are needed to find the three quantities A™ = q„i{<J^) / o'^ {m = 0,1,2), as well as the parameter 
K {rj, t, a). We stress that the almost fully analytical determination of these unknown parameters was lacking in Ref. 
[32] and represents an important part of the present work. Our detailed analysis is given in Appendix B, and we 
quote here the main results. 

i) For Aq, the same PY equation found by Baxter for isotropic SHS [ill. Il2| 

I2r,t Ag - (^1 + i^f^ Ao + 2/f (r/)i = 0. (5.1) 

Only the smaller of the two real solutions (when they exist) is physically significant [III, [I^, and reads 






12lLt+Jli + j2iLty-48vy^Hv)t' 



(5.2) 
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1277,„i A^ - ( 1 + -3^^t ) A,„ + /i^^(77„Oi = -V (m = 1, 2). (5.3) 



ii) For Ai and A2, two other quadratic equations, i.e. 



,2 / 1 , 1277,„ > HS/ 



iii) The fourth equation is the foUowing hnear relationship between Ai and A2 

12?72 A2 12?/i Ai _ r/2 (4 - 7^2) ??! (4 - jyi) 



i-'72 1-771 (1-772) [l-ViY 



(5.4) 



which stems from the condition a2 = ai. 
The analysis of Appendix B gives 



with 



A2 (771,772,^,0) = Ai (7/2, 771, t, a) (5.5) 



A,„=A + A^ (™=1,2) (5.6) 

A :. 1 + 1 (^]^ + -^]^] = i + -^(1 + 4^) f5 7) 

3^4^1-7/1 1-772/ 3^4(l + a;)(l-2a;) ^^^ 

Aox _ V^ „^ex AGX _ ^1 H^ex /r Q\ 

' ~ ^{l-V2) ° ' ' "4(1-771)^0 ' ^^■^> 

where we have introduced 771 = —a;, 772 = 2a; { x = Krj ), and WJ'^ is defined in Appendix B. All these quantites are 
analytic functions of a; = /\77. Thus, to complete the solution, we need an equation for K^ which can be written as 

(5.9) 



K = at /C, 


with /C. f<"> , 

^(771, 772, i) 


= ^(Ai+A2)-3|i>; 


'l277™A^-i^^=^ + /7H^(7;„0" 



^Y (5.10) 

and lim,;^o Z{i]i, 772, t) = 1. Insertion of found expressions for Ai, A2 and A',og (sec Appendix B) into Eq. (|5.9p yields 
a single equation for K that we have solved numerically, although some further analytic simplifications are probably 
possible. 

Our solution is then almost fully analytical, as only the final equation for K is left to be solved numerically. 

B. Approximate expressions 

For practical use we next derive very accurate analytical approximations to K^ Ai and A2, which provide an useful 
tool for fully analytical calculations. Since in all cases of our interest we always find x = Kr] <C 1, a serie expansion 
leads to: 

W^^ = ^{l + 5x)t + 0{x^), (5.11) 

and, consequently, 



Similarly we can expand Z in Eq. (|5.10p as 



Z{x,t) = 1 + zi{t)x + Z2{t)x'^ + O (x^) 



with 



zi(t) = |(3 + llt) 

Z2{t) ^ I {15 + 6lt ~ Af) . 

Insertion of this result into Eq. (|5.9p yields a cubic equation for K, 

Z2{t)TfK'' + Zi{t)ijK^ +K-at yPY(a) = 0, 

which, again with the help of Eq. (|5.9p . is equivalent to a cubic equation for Z 

Z^ -Z' + z,{t) [at yr{<^)v\ Z + Z2{t) [at yr iM ' = «■ 
The physically acceptable solution then reads 



Z(77,t) = i [i+\Jb+^B^^C^+ '\Jb - v/S2 - C3 j , 



where 



6=1- 
C = l 



|zi(i) 
3zi(i) 



at yl^{a)-q_ 
at yl^{a)r] 



^:Z2{t)[atyr{M- 
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(5.13) 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



In conclusion, our approximate analytic solution for K, Ai and A2 includes three simple steps: i) calculate K by using 
Eqs. dEH), ((5TT6l) - ((CT7l) . ((511)1 : ii) evaluate x = i^T?; iii) solve for Ai and A2 by means of Eqs. ([5?7|) and ([5l^ . 



C. Numerical comparison 



In order to assess the precision of previous approximations, we have calculated K^ Ai and A2 by two methods: i) 
solving numerically Eqs. (|B8p . and ii) using our analytic approximations. After fixing a — 1/2, we have increased 
the adhesion strength (or decreased the temperature) from t = (HS limit) up to f = 0.8, for some representative 
values of the volume fraction (77 = 0.01, 0.1, 0.2 and 0.4). The maximum value of t corresponds to r = l/(12t) ~ 0.1, 
which lies close to the critical temperature of the isotropic SHS fluid. On the other hand, 77 = 0.01 has been chosen 
to illustrate the fact that, as 77 -^ 0, the parameter K tends to at. The linear dependence of K on t in this case is 
clearly visible in the top panel of Figure 2. 

In Figures 2 and 3 the exact and approximate results for K , Ai and A2 are compared. The agreement is excellent: 
at ?7 = 0.1, 0.2 and 0.4, the relative error on K does not exceed 0.1%, 0.4% and 1%, respectively, while the maximum 
of the absolute relative errors on Ai and A2 always remain less than 0.05, 0.2 and 0.6 % in the three above-mentioned 
cases. It is worth noting that, as rj increases, the variations of Ai and A2 are always relatively small; on the contrary, 
K experiences a marked change, with a progressive lowering of the relevant curve. 



VI. SOME ILLUSTRATIVE RESULTS ON THE LOCAL ORIENTATIONAL STRUCTURE 

Armed with the knowledge of the analytic expression for the q^s a rapid numerical calculation of the three harmonic 
coefficients {/iq, /ia, ho} appearing in 
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PY-OL 



(1, 2) = 1 + ha{r) + /iA(r)A(l, 2) + hD{r)D{l, 2). 



(6.1) 



can be easily obtained as follows. From the second Baxter IE (|4.8p . one can generate h{r) directly from q{r)^ avoiding 
the passage through c{r). From {90191192} one first obtains {/ig, /ii, /i2}, by applying a slight extension of Perram's 
numerical method |42| and then derive {/iq, ^a, /in}, according to the above-mentioned recipes. 
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The main aim of the present paper was to present the necessary mathematical machinery to investigate thermo- 
physical properties. We now illustrate the interest of the model by reporting some preliminary numerical results on 
the orientational dependence of ^^^"'-"^(l, 2) - i.e. on the local orientational structure - as a consequence of the 
anisotropic adhesion. A more detailed analysis will be reported in a forthcoming paper. 

Consider the configuration depicted in Figure 4. Let a generic particle 1 be fixed at a position ri in the fluid with 
orientation ui , and consider another particle 2 located along the straigth half-line which originates from the center 
of 1 and with direction Ui . This second particle has then a fixed distance r from 1, but can assume all possible 
orientations U2 , which - by axial symmetry - can be described by a single polar angle 9 = 02 (i.e., the angle between 
Ui and U2 ) with respect to the intermolecular reference frame. Within this geometry, we have {9i,(pi) = (0,0) and 
</j2 = 0, obtaining A(l, 2) = cosO, D{1,2) = 2cos0. Consequently, 5(1,2) = g{r,9i,ipi, 62,^2) reduces to 

g{r,6)^gQ{r) + [h^{r) + 2hD{r)]cose, (6.2) 

where 9 = 62, and goir) = 1 + hoir) is the radial distribution function of the reference isotropic SHS fluid. 

Clearly, g{r, 9) is proportional to the probability of finding, at a distance r from a given molecule 1, a molecule 2 
having a relative orientation 9. We consider the three most significant values of this angle: i) 9 = Q, which corresponds 
to the 'parallel' configuration of ui and U2; ii) 9 = 7r/2, for the 'orthogonal' configuration; and 9 ^ tt, for the two 
'antiparallel' (head-to-head and tail-to-tail) configurations. From Eq. (|6.2p it follows that 

gP-(r) = g{r, 0) = go{r) + [h^{r) + 2hD{r)] , 

gortho(^)^^(^^^/2) = go(r), (6.3) 

^a„tipar(^) = g{r, tt) = go{r) - [/iA(r) + 2hD{r)] . 

Note that g°'^'^°(r) coincides with the isotropic result go{r). 

In Figures 5 we depict the above sections through the three-dimensional surface corresponding to g{r^ 9), i.e., g^^^{r), 
gOTtho^^^ and ^antipar^j,^^ £qj. ^ _ q g -^^jth t — 0.2 and t ~ 0.6, respectively, at the highest asymmetry value admissible 
in the present model, i.e. a = 1/2. The most significant features from these plots are: i) g™*'P^'((7+) > gP'^''(cr+); ii) 
for r > 2(7 (^antipar^j,-| ^ gpni ^^^ ^ go{r), i.e., the anisotropic adhesion seems to affect only the first coordination layer, 
a < r < 2(7, around each particle. 

The interpretation of these results is the following. In view of i) we see that the parallel configuration is less probable 
than the antiparallel one at contact. Such a finding, together with ii), means that chain formation characteristic of 
polymerization is inhibited by the short-ranged anisotropic adhesion exploited here. This strictly contrasts with the 
case of long-ranged DHS fiuids, where it is believed [7|, |8| that chaining phenomena might preempt the gas-liquid 
transition. This specific feature of the present model is extremely interesting and we plan a throughout investigation 
on this topic in a future publication. 



VII. PHASE STABILITY 

In view of previous findings, a very natural question is whether the addition of our anisotropic sticky term to the 
potential changes phase stability and phase transition curves with respect to the corresponding isotropic case. We 
believe the answer to be positive. This is stron gly suggest ed by results obtained for similar anisotropic models, such 
as hard spheres with 'sticky points' [ll,[ll,|ll,lil,[llli,[li,0,[ll| or 'sticky patches' 0, [H, IH, |2l, [H, HI, Hi . 

We now briefly comment on this issue. Within our formalism, this problem of stability can be conveniently analyzed 
using standard formalism devised to this aim [iji Ea, H, |43| . 

We start from the stability condition with respect to small but arbitrary fluctuations of the one-particle density 
p(l) from the equilibrium configuration, denoted as 'eq' |45l. |46|. |47| . 



d(l) / d{2) 



<5(1,2) 



P(l) 



c(l. 



Sp{l)Sp{2) > 0. (7.1) 



eq 



Here d(i) stands for dvi dQ,i, i = 1,2, and we assume the equilibrium one-particle density to be p/Air [45ll46l. l47j . 
We expand the fiuctuations both in Fourier modes and in spherical harmonics |44l | 

/,, -l-oo +1 

7:rT3 e*-""^' T.T.^P (k) Yim (n,) . (7.2) 



/=0 m=-l 



Using the orthogonality relation [44 
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standard manipulations [47| show that condition (|7.ip can be recast into the form 



(7.3) 



dk 



+00 +h +I2 

51 XI XI / 77-73 '^^'i'"! (^)'^^;*i™l W^'i™i'2rn2 (k) >0 

ll.l2—0mi— — li 7712 — — ^2 ^ "^ 



where the matrix elements Ai^rnii27n2 (k) are given by 



(7.4) 



An 



X f dr&^''c{r,ni,n2). 



dill / dr!2l^,^, (f^OF;,™, (J72) 



(7.5) 



The problem of the stability has been reported to the character of the eigenvalues of matrix (|7.5p . This turns out 
to be particularly simple in our case. Using the results (jA3p it is easy to see that 



(7.6) 



/'drc'^-"c(r,r!i,r!2) = £o (fc) + Sa (fc) A (f^i, f^a) + cjj (fc) i? (f^i, ^2, ^ifc) 
Insertion of Eg. ([7^ into Eg. (173)1 leads to 



47r 

■^limil2'm,2 V^J V -^J ^lilri^nii. — 7n2 

P 



(7.7) 



Co (fc) liZ-.hrm + ^A (fc) ^itnii^ms + ^^ (^) ^imii^ma 



where we have introduced the following integrals, which can be evaluated in the intcrmolecular frame, using standard 
properties of the spherical harmonics [44| 



/, 



/, 



(0) 
Iimil2m2 

(A) 
Iiniil2va2 



= / dOi / dil2 Yi^rni {fll)Yi^m2 {^2) A (ill, 1^2) = -;:TrSl^iSi^^iS,nioS: 

4mii2m2(cos^) = JdniJdn2Yi,^,i^i)Yi,^A^2)D{ni,n2,nk) 



1712,0 



= oTTfJzi^iiJi^, 1(5^1, o(5„2,o 2 P2 (cos^) 

and where ^2(2;) = (3x^ — l)/2 is the second Legendre polynomial. 
Hence, the matrix (|7.5p is diagonal and the relevant terms are 



(7.8) 



^0000 (k) = Air 



Co (fc) 

P 



whose positiveness is recognized as the isotropic stability condition, and 

1 1 



^1010 (k) = 47r - - - [£a (fc) + 2P2 {cos 0) CD (fc)] 
.P 3 



(7.9) 



(7.10) 



All remaining diagonal terms have the form Aiqiq — An / p > 0. 
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In order to test for possible angular instabilities, we consider the limit fc -^ of Eq. (|7.10|) namely 



iioio(O) = — {l-^[cA (0) + 2P2 {cose) com}- (7.11) 



This can be quickly computed with the aid of Eqs. (tSUOl), (|^ - the fact that cd{0) = c%{0) and the identity (IJTTS)) . 
We find 

Aioio(O) - —al (7.12) 

which is independent of the angle 9. This value is found to be always positive as ai > (see FiglH]). Within 
this first-order approximation, therefore the only instability in the system stems from the isotropic compressibility. 
The reason for this can be clearly traced back to the first-order approximation to the angular dependence of the 
correlation functions. If quadratic terms in A and D were included into the series expansion for correlation functions, 
the particular combination leading to a cancellation of the angular dependence in the stability matrix Ai-^„iii2m2 (0) 
would not occur, leading to a different result. 

This fact is consistent with the more general statement that, in any approximate theory, thermodynamics usually 
requires a higher degree of theoretical accuracy than the one sufficient for obtaining significant structural data. 
Conceptually, the need of distinguishing structural results from thermodynamical ones is rather common. For instance, 
in statistical mechanics of liquids it is known that approximating the model potential only with its repulsive part 
(for instance, the hard sphere term) can account for all essential features of the structure, but yields unsatisfactory 
thermodynamics. On the other hand, the present paper refers to a simplified statistical-mechanical tool, i.e. the OZ 
equation within our PY-OL closure, which has been explicitly selected to allow an analytical solution. Our results 
however indicate that the first-order expansion used in the PY-OL closure can give reasonable information about 
structure, but not on thermodynamics, where a higher level of sophistication is required. 

VIII. CONCLUDING REMARKS 

In this paper we have discussed an anisotropic variation of the original Baxter model of hard spheres with surface 
adhesion. In addition to the HS potential, molecules of the fluid interact via an isotropic sticky attraction plus an 
additional anisotropic sticky correction, whose strength depends on the orientations of the particles in dipolar way. 
By varying the value of a parameter a, the anisotropy degree can be changed. Conscquenly, the strength of the total 
sticky potential can vary from twice the isotropic one down to the limit of no adhesion (HS limit). These particles may 
be regarded as having two non-uniforrn, hemispherical, 'dipolar- like patches', thus providing a link with uniformly 
adhesive patches [H, M, HI, Hi, [H, [H, [S^ . 

We have obtained a full analytic solution of the molecular OZ equation, within the PY-OL approximation, by using 
Wertheim's technique [3|. Our PY-OL approximation should be tested against exact computer simulations, in order 
to assess its reliability. Nevertheless, we may reasonably expect the results to be reliable even at experimentally 
significant densities, notwithstanding the truncation of the higher-orders terms in the angular expansion. Only one 
equation, for the parameter K, has to be solved numerically. In additon, we have provided analytic approximations 
to K, Al and A2 so accurate that, in practice, the whole solution can really be regarded as fully analytical. From this 
point of view, the present paper complements the above-mentioned previous work by Blum et al. [3^ . 

We have also seen that thermophysical properties require a more detailed treatment of the angular part than 
the PY-OL closure. Nonetheless, even within the PY-OL oversimplified framework, our findings are suggestive of a 
dependence of the fluid-fluid coexistence line on anisotropy. 

Our analysis envisions a number of interesting perspectives, already hinted by the preliminary numerical results 
reported here. It would be very interesting to compare the structural and thermodynamical properties of this model 
with those stemming from truly dipolar hard spheres [4^, |4y, |43] . The possibility of local orientational ordering can be 
assessed by computing the pair correlation function g{l, 2) for the most significant interparticle orientations. We have 
shown that this task can be easily performed within our scheme. This should provide important information about 
possible chain formation and its subtle interplay with the location of the fluid-fluid transition line. The latter bears a 
particular interest in view of the fact that computer simulations on DHS are notoriously difficult and their predictions 
regarding the location of such a transition line have proven so far unconclusive [ij] . The long-range nature of DHS 
interactions may in fact promote polymerization preempting the usual liquid-gas transition [8|. Our preliminary 
results on the present model strongly suggest that this is not the case for sufficiently short-ranged interactions, thus 
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allowing the location of such a transition line to be studied as a function of the anisotropy degree of the model. Our 
sticky interactions have only attractive adhesion, the only repulsive part being that pertinent to hard spheres, whereas 
the DHS potential is both attractive and repulsive, depending on the orientations. 

Finally, information about the structural ordering in the present model would neatly complement those obtained 
by us in a recent parallel study on a SHS fluid with one or two uniform circular patches [10| . Work along this line is 
in progress and will be reported elsewhere. 
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APPENDIX A: EXTENSION OF WERTHEIM'S APPROACH 

The Fourier transform of the excess PY-OL equation, Eq. (|3.11[) . reads 



/lex(k,r2i,f72) =2cx(k,r2i,f72)+P ( 2cx(k,r2 1,173) /lox(k,f73,^2) ) (Al) 

(the superscripts have been omitted for simplicity). In order to evaluate the angular average, we first need the FT of 
c and h. The FT integral (|3.2[) may be rewritten as 

/ drr^ dilr exp(ik-r) [•••]=/ dr r^ d(j> d{cos9) e''^''''"^'^ [■■■]. 

Jo J Jo Jo J-1 

Let us now apply this operator to i^ox(l, 2) (F = c, h), expressed as 

Fex(r,0i,r!2)-i^A(r) A(r!i, f^a) + i^cW D(17i, ^82, ^i.), (A2) 

and first perform the angular integration J di^r, recalling that [3| 

/ drir exp(ik • r) 1 = 47r jo{kr) 1 
J drir exp(ik • r) A(f7i, 172) = 47r jo{kr) A(fii, O2) (A3) 

JdQr exp(ik-r) /^(rii, 172, ^r) = -47r J2(fc?') ^(^1.^2,^/=), 

where joix) = x^^ sinx and J2ix) = 3x^^ sin a; — 3x~^ cosx — .70(2;) are Bessel functions, and 

£)(17i, 172, rjfc) = 3(ui • k)(u2 • k) - ui • U2 = -Dfc(l, 2), 

with k = k/fc. We get 

Pcx(k,17i,172)-^A(fc) A(17i,172)+^D(fc) i?(17i,172,17fc), 

where F^ik) is the usual FT of the spherically symmetric function FA(r): F,,,{k) ~ 47r Jj^ dx x"^ jo(fcx)F..(x-). On 
the other hand, F ij{k) ~ — 47r J„ dx x^ J2(kx)F£){x), which is the Hankel transform of Fjjir), may conveniently be 
considered as the FT of a 'modified' function F'l){r), i.e. Foik) ~ F^{k). Taking the inverse FT of Foik) yields 

F^D {r) = ^J°°dke JO {kr)FD (k) ^ Foir) ~ 3 J^ ^^ dx, (A4) 

with the help of the identity 

W {x — r) S {x — r) 



dk k jo{kr) J2{kx) = - 
^ 
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In conclusion, the FT of Fcx(l, 2) reads 

Pcx(k,r!i,f]2) = ^A(fc) Aini,n2) + FUk) Dini,n2,nk), (A5) 

with F standing for h or c. 

Let us now define the angular convolution of two functions as 

AoB^BoA={A{ni,n3) B{n3,n2)),,^. 

Wertheim [3| demonstrated that the rotational invariants 1, A, and D form a closed group under angular convolution; 
that is, the angular convolution of any two members of this set yields only a function in the same set, or zero, according 
to the following multiplication table 



A D. 



110 

A A/3 Dk/S 

Dk Dk/3 (Dfc + 2A)/3 

TABLE I: Angular convolutions of the basis functions 1, A and Dk- 



Substituting the expressions for Cex and /lex given by Eq. (|A5[) into the angular average Ccx ° ^cx = 
Cox(k, 1^3, 5I2) /iox(k, fls, 1^2) ) , with the help of Table I we obtain 

Ccx o Kx = ca/ia o ^ + c/^h%-Dk 

+7Pji^'^Dk + ^dKI (2A + Dk) . 

Inserting this result into Eq. (jAip and equating the coefficients of A, and D separately, one finds that the k-space 
excess PY-OL equation splits into two coupled integral equations, i.e.. 

Ha -CA = hp ( caIt-a + 2 c^^/i^) 



^H \ '-A"'A -T- ^ ^D"'D 

I 

3 

Coming back to the r-space, one gets the following equations 



h'h-^D = \p CAh%+7^\,hA+^hhl 



I /ia(?-) = CA(r-) + ip ( CA^/iA + 2 c?)*/!?, ) ,^^. 

\h%{r)=clir) + y {cA*h%+cl*hA + c»D*h%). 

In particular, since h£i{r) = for < r < cr, Eq. (|3.27p yields h'-j^{r) ~ — 3A' for < r < cr, with K being a 
dimensionless parameter defined by 

K^r^-^d.. (A7) 



The exact core conditions for Eqs. (|A6|) are 

Y'l^^ ] forO<r<a. (A8) 

Now, in the PY-OL closure for the DCFs, Eqs. (|3.18p . the closure for coir) must be replaced with that corresponding 
to c'jj{r) (for simplicity, here and in the following we omit the superscript PY-OL). In order to derive this, let us start 
from Coir) = CD,rcg{r) + Ad <tS (r - ct) , where CD,rcg{r) = f^^{r) y^ {r) = for r><7. Then Eq. ([3?27|) yields 

cl{r) = CD{r) - 3 f ££££iM dx - SAd 9{a - r), 
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since J^ S {x ~ a) x~^ dx ~ a~^6{a — r) [S^- So we get 

Cc(r) ==C£,(r) for r > a, (A9) 

and the required new closures are 

c%{r) = Ad a6 {r - a) f " ^ ' 

In order to decouple the two integral equations for A- and Z?-cocf5cients, we then introduce two new unknown 
functions, which are linear combinations of the previous ones. Defining 

-Fnew = AiFa + X2F^ {F = C, h) 

and using Eqs. @ leads to 



AiCA^A + A2 (cAho + ^d^a) + (2Ai + A2) '^D^' 



1 



3' 

Requiring the second member of this equation to be proportional to pcncw^now ^ that is, equal to 
Cp (Aica + A2c2j) I Aift-A + ^2^% ) , with C being the proportionality constant -, yields the following conditions 

3A2 = CX1X2 

i(2Ai+A2) = /:A2 . 
An infinite number of solutions are possible, and correspond to 

(Ai,A2) = ^(l,-l), and (Ai, A2) = ^ (1, 2) , 

since there is no need for the proportionality constant to have the same value in the two cases, i.e. ^2 can differ from 
^1. As a consequence, we can write the two new hacwir) as 

f h,{r)^i3C^)-'[hAir)-hUr)] ..... 

\h2{r)^{3C2r'[hA{r)+2hl{r)], ^ ' 

while similar expressions hold for ci and C2. From Eqs. (|A8[) it follows that: hi{r) = K/Ci and h-zir) ~ —2K/C2 
for < r < (J. 

In Ref. I Wcrthcim chose £1 = —K and £2 = ^K Q, which leads to 

for < r < cr (A13) 

(in Ref. I, Fi and F2 were denoted as F_ and F+, respectively). Clearly, Wertheim's choice has the advantage of 
providing, for all the three 'hypothetical' fluids, core conditions of the typical HS form: hm{r) = — 1 for < r < cr 
(to ~ 0, 1, 2). The cost to pay is the introduction of 'modified densities' for the 'auxiliary' fluids 1 and 2 (the negative 
sign of pi poses no special difficulty). 

On the other hand, it is would be equally proposablc the choice £1 ~ £2 = 1, which leads to 

[F,{r)^\ [FA(r)-F^(r)] 
\F2[r) = \ [FA(r)+2i^0(r)] 

\P2^P, \h2{r)=-2K 

The advantage of this second possibility would be that all the three 'hypothetical' fluids have the same real density, 
while the cost is represented by the less usual core conditions, which however pose no particular difficulty. 
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APPENDIX B: EQUATIONS FOR THE UNKNOWN PARAMETERS 

Three quadratic equations for the A„i 's (m = 0,1,2) can be obtamed from Eqs. (|4.3p - (|4.4p . after deriving from 
Eq. (j4.17p the following expressions for the PY-OL contact values 



2 



where 



hm,rcg (ct+) = h^^ (r;„0 - -i^^^ A,„ + 12?7,„ A: 

yHS (^) = (1 + 1^) (1 _ :,)-2 . 

Substituting Eq. (jBip into the expressions for A^ given by Eqs. (|4.3p , we get: 
i) for Aq, the same PY equation found by Baxter for isotropic SHS [lll.[ia| 

12vtAl-^l + ^^t^Ao + y^Hv)t^O. 

Only the smaller of the two real solutions (when they exist) is physically significant [III [121 1 and reads 



(Bl) 



(B2) 



(B3) 



An 



ii) For Ai and A2, the equations 



1277,„i Al - 1 



■Ja^'iv) t 






12??,n 
1 -?7n 



1277 



l + -f^ 



487? 2/HS(,^) i2 



t A„ + K^{ri,n)t = -V {m = 1, 2). 



(B4) 



(B5) 



It is remarkable that the right-hand member of these equations does not depend on the index m. This fact means 
that A2 obeys exactly the same equation as Ai, but with 772 replacing 771; as will be confirmed later, such a property 
implies that, if one writes Ai = Ai (771,772, t, a), then A2 must have the same functional form with 772 interchanged 
with 7^1 , i.e. A2 (771 , 772 , i, a) = Ai (7^2 , 771 , t, a) . 

Now the system of equations for Ai, A2 and K must be completed by a further relationship, which can be obtained 
from the sum rule, Eq. (|3.46p . Taking into account that c% = K (2c2 + ci), and multiplying Eq. p.46p by Airp yields 



47rp2 / C2{x) X dx ^ 47rpi / ci{x) x dx, 
'0 Jo 



(B6) 



On the other hand, puttmg k = into Eq. (|4.7p gives 

poo 

1 - PmCm (fc = 0) = 1 - 4:TTpm / Cm{r) r^ dr = Ql^{k = 0) = a^ 

since Qm{k = 0) = a„i (as shown by the first of Eqs. (|4.10p ). Then Eq. (jB6P becomes a: 
two equations: 02 = ai, and 02 = —ai. From the expression for am, one can easily realize that the second equation 
does not satisfy the i ^ limit, whereas the first one. 02 = ai (or, equivalently, a"'°^^^(772, A2) = a"'°^^^(77i, Ai) ), 
leads to the following linear relationship between Ai and A2 



af, which splits into 



12772 A2 12771 Ai 



HS/ 



m 



Vi 



a^^m) a"^(r;i) 



(B7) 



Note that the two Eqs. ([BS]) are coupled (since Krcg/K = 1 - Ad/K = 1 - (2A2 + Ai) ), but with the help of Eq. 
B7p they could be easily decoupled. However, since the right-hand members of Eqs. (jBSp coincide, we can get a new 
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relationship by equating their first members, and exploiting Eq. (jB7[) . So we arrive at the following equations for the 
three unknowns Ai, A2 and K: 



12ri2t KI-K2+ 6HS(^2)i = 1277ii A? - Ai + b'^'^{rii)t 

12r)2 A2 _ 12171 Ai _ r)2 (4-1)2) _ i)i(4-t;i) 
l-r;2 1-1)1 ~ (1-1)2)^ (l-i)i)'' 



1277ii A? - (1 + i22^t) Ai + h^''{m)t - -^- 
The first two equations form a closed system for Ai and A2. The second one suggests that we can assume 



(B8) 



1277„i A„i _ 77„i (4-77„,) 



1 - n,n {i-nn,y 



W, 



or, equivalently, 



A„, = 



1 , "Hn 



1 -77n 



W (m = l,2). 



(B9) 



3 4(l-7/,„) 12r;„ 

where W = M^ (771, 772, i) is an unknown function, which must be proportional to 771772- In fact, Eqs. (jBSp require that 

1 



lim Ai ~ lim A2 — , 

r;^0 1)^0 3 



(BIO) 



since, from Eq. (|4.4p . one has lim^^o ^ ~ \ ^ hm,,_,o /"^Ticg = 0). If Ai and A2 in the first of Eqs. (|B8p are replaced 
with the new expressions (|B9[) . then one gets a quadratic equation for W: 



with 



where we have put 771 



with 



(1 - 771772) iW^2_(^_ 2771 7,2t)W^ 



3771772 



(1 - 7;i) (1 - 772) 



M = 



M = 1 
= 1 



1 + 2771+772)^5771772 1 

Tj TTJ ^ o 1 ~ '7l) 1 - m) 

(1 - 7;i) (1 - 772) 3 

1 + 2a; + 10x2 1 , ^^_, ^ J ^ 

- -(l + a;)(l-2a:) i. 



_(l + a;)(l-2x) 3 
-x, 772 = 2a; ( a' = ^^77 ). The acceptable solution is 



2 (1- 7717^2) t 
37/1772 



Wn 



(1 - 7/1) (1 - 7/2) 



M 



Wn = - 



6x 



(l + x)(l-2a;) 



Wo 



M 



i(l-277i7/2t)(l + VP) ^{l + AxH)(l + VV 



(Bll) 



(B12) 



(B13) 



(B14) 



V = 1 



127717/2 (1 - 771772) 



(1 - r/i)(l- 7/2) (1- 27/1 772i)' 



Mt 



24a;2 (l + 2a;2) 



(l + a;)(l-2a;)(l+4a;2t)' 



Mt. 



Note that hm,,,^o Wq = lim,,^o il/ = 1 + (2/3) i. 



(B15) 



The functions W, Wq, D and M are symmetrical with respect to the exchange of 771 and 7/2; in particular, 
1^(772,771,^) = 14^(771,772, i), and this property implies that 



^2iVl,V2,t) = Ai (772, 7/1, t). 



(B16) 
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confirming our previous guess. 
Moreover, if we put 



then 



with 



A: 



1 1 



m 



M^o = 1 + Wo^^ 



A — A + A"'' 



V2 



x{l + Ax) 



3 4Vl-'/i l-??2/ 3 4(l + a;)(l-2x)' 



(B17) 
(B18) 

(B19) 



Af = 



V2 



4(1-772) 



H^o^ 



^r = TTT^^W^""- 



4(l-r?i) 



(B20) 



Here, both A and W§^ are symmetric with respect to 771 and 772: whereas A^ represents the asymmetric part of A^- 
Note that the knowledge of Ai and A2 allows to calculate Aa and Ad immediately. In fact, Eqs. (|3.4ip lead to 



Aa = 2K (A2 - Ai) = -K 2(1+4^1-2.) W§- 

Ad = /^(2A2 + Ai) ^k[i+ ,^,J{^,_,,^ [1 + 2a; (2 + M^o-)]} 



(B21) 



Now we must find an equation for K. We can regard the third of Eqs. (jB8[) as the required relationship. However, 
in order to derive a more symmetric expression, we prefer to start from Eqs. (jB5[) . rewritten as 



and we get 



1277it A? - (1 + ^A Ai + h^^Mt + %t + if yVi'^) = 
12772* Ai - (1 + ii^t) A2 + h^^^,)t + %t + if yna) = 0, 



z{m,V2,t)' 



K = at K., with K, = 



(B22) 



(B23) 



Z=^(Ai+A2)-3 i^ 



1277™ K 



1277,„A,,i 

1 - T^m 



+ /ir(^/™) 



Kr, 



K 



(B24) 



and lim^^o Z{rii, 772, i) — 1. Replacing the found expressions for Ai, A2 and A^i into Eq. (|B23p yields an equation for 
K that wc have solved numerically, although some further analytic simplifications are probably possible. 
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a) Isotropic adhesion 




b) Anisotropic adiiesion: 
parallel configuration 





c) Anisotropic adhesion: 
antiparallel configurations 



FIG. 1: (Color online) Illustration of the dipolar-like adhesion. In the top panel a) the adhesion is isotropic, with e(l, 2) — 1. 
In the two other cases the adhesion is anisotropic and: i) stronger and maximum in the head-to-tail parallel configuration 
b), where e(l,2) = 1 + 2q; ii) weaker and minimum in the two antiparallel configurations c) ( head-to-head and tail-to-tail 
orientations, both with e(l, 2) = 1 — 2a). 
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FIG. 2: (Color online) Comparison between exact numerical and approximate analytical results for the parameters K, Ai, and 
A2 as a function of t, for anisotropy degree a — 1/2 and two values of the packing fraction: 77 = 0.01 (top panel) and 77 = 0.1 
(bottom panel). 
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FIG. 4; (Color online) Illustration of the simple configuration discussed in the text and chosen to define some radial sections 
through the multidimensional plot of g{l,2). 
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FIG. 5: (Color online) Sections through (;(1,2), with particles 1 and 2 in the configuration shown by the previous Figure, 
calculated as a function of r for fixed relative orientations: = (parallel configuration), Q = 7r/2 (orthogonal configuration), 
and 9 = TT (antiparallel configurations). 
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FIG. 6: (Color online) Evaluation of quantities ao (top panel) and ai — a2 (bottom panel) as a function of t for various packing 
fractions ranging from 77 = 0.01 to r; = 0.4. These are computed from Ea. (|4.12p with m = 0, 1. Note that for both 77 = 0.1 and 
rj = 0.2, ao — corresponds to the onset of isotropic instability. 



